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Abstract 

The review of the following results of the Refs, [lj - [5] is presented: For mixed state 
light of TV- mode electromagnetic field described by Wigner function which has generic Gaus- 
sian form the photon distribution function is obtained and expressed expliciltly in terms of 
Hermite polynomials of 2TV-variables.The momenta of this distribution are calculated and ex- 
pressed as functions of matrix invariants of the dispersion matrix.The role of new uncertainty 
relation depending on photon state mixing parameter is elucidated.New sum rules for Her- 
mite polynomials of several variables are found.The photon statistics of polymode even and 
odd coherent light and squeezed polymode Schrodinger cat light is qiven explicitly.Photon 
distribution for polymode squeezed number states expressed in terms of multivariable Her- 
mite polynomials is discussed. 


1 Introduction 

In the Ref.fl] it was shown that the matrix elements of density matrix in number state basis for 
polymode oscillator are expressed in terms of Hermite polynomials of several variables for the 
density operator in the canonically transformed thermal state of the oscillator.In the recent works 
[2], [3] the photon distribution function for the generic Gaussian light described by the Wigner 
function which is the most generic Gaussian in quadrature phase space was found and expressed in 
terms of Hermite polynomials of 2 variables for one-mode case [2] and 2N variables for polymode 
case [3]. The physical meaning of mixed Gaussian state of the light may be understood if one takes 
into account that the pure multimode Gaussian state corresponds to the generalised correlated 
state introduced by Sudarshan [6] who related those states to the symplectic dynamical group.The 
mixed Gaussian states studied above may be considered as the mixture of generalised correlated 
states plus thermal noise acting on each mode which has its own temperature. The photon 
distribution function for even and odd coherent states [7] or Schrodinger cat states [8] subject 
to squeezing both in one- mode and polymode cases has been found in Ref. [4]. The polymode 
Schrodinger cat states and photon distributions for light in these states were introduced in [5].The 
aim of this work is to give a review of the photon distribution functions and related sum rules 
for one-mode and polymode Gaussian light and for the even and odd coherent states light using 
the results of [2] - [5]. It should be emphasized that the squeezed light is worth to be used in 
interferometric gravitational antennas [9] and the even and odd coherent light may play alternative 
role in gravitational wave experiment [10], [11]. 
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2 Photon Distribution for Polymode Gaussian Light 


The mixed squeezed state of the iV-mode light with a Gaussian density operator g is described 
by the Wigner function (see, for example, [12]) 

W’(p.q) = (detM)-lexp [-|(Q- < Q >)M" 1 (Q- < Q >)] , (!) 

where 2 J /V-dimensional vector Q = (p,q) consists of N components and N components 

q u ..., q N . 2 N parameters < Pi > and < « >, * = 1, 2, . . . , N, combined into vector < Q>, are the 
average values of the quadratures. A real symmetric quadrature dispersion matrix M has 2N +N 
parameters 


M a 0 = ~ (QaQp + QpQa) ~ (Qa) {Qp) » <*,/? — 1 , 2 ,..., 2JV. (2) 

They obey certain constraints, which are nothing but the generalized uncertainty relations [12]. 
The photon distribution function in this state has the form [3],[1] 


V Q = Vo 


n! 


(3) 


where vector n consists of N nonnegative integers: n = (ni,n 2 ,...,njv)- The function Hnn (y) 
is the Hermite polynomial of 2 N variables. We introduced also notations 

n! = n!!n 2 !...n^!. (4) 


The symmetric 
relations 


2/V-dimensional matrix R and the 2iV-dimensional vector y are given by the 
R = U f (1 2N - 2M) (lajv + 2M) _1 U% (5) 


y = 2 U*(I 2 n - 2M) _1 < Q >, (6) 

where the 2iV-dimensional unitary matrix 

tt — — ( ~^ N 

y/2\ In In ) 

is introduced. The matrices Ijv and I 2 n are identity matrices of corresponding dimensions. 

The probability to have no photons has the form 

Vo = [det * exp [- < Q > (2M + I 2 n ) -1 < Q >] - ( 7 ) 

It may be shown [3], [4] that the multivariable Hermite polynomial is the even function if the 
sum of its indices is the even number and the polynomial with the odd sum of indices is the odd 
function. Due to this pairity property of the polydimensional Hermite polynomials the ’’diagonal” 
multivariable Hermite polynomial is the even function since the sum of its indices is always even 
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number.Consequently the above photon distribution function is the even function. Thus the 
photon distribution function for generic mixed Gaussian light found in[3] is expressed in ter ms of 
multivariable Hermite polynomials and it depends on the quadrature means and dispersions. The 

photon number means and dipersion matrix corresponding to the found distribution (3) are of the 
form 

< n : >= - 1) + i(< pj > 2 + < q J >), 

an i n j ~ 2 ~ 2 dj - -)+ < Qj > Mj < Q j >, 

where Tj and dj are the trace and the determinant of the 2x2- matrix Ad j, describing only j-th 
mode, and the 2- vector Q j has the components ( Pj,qj )• 


3 Pure Polymode States 


The photon distribution for polymode squeezed correlated state may be expressed in terms of 
symplectic transform parameters relating boson operators as follows 



( 8 ) 


where ft is a symplectic 2 A r x2. , V- matrix consisting of four TV-dimensional complex square blocks, 
and d is a complex IV-vector.Then we have for photon distribution in squeezed correlated polymode 
state |/? > labeled by the complex number vector with IV-components 

p„ = ^)| W K-'.)(,-.[ /3 _ d ])| I > Vo(fi) — \Fo(fl)\ 2 exp —\/3\ 2 t (9) 

where 

W) = (Jet cr* exp [j/VC ‘/3 + w - iT‘d) + itVC'd - i|d| J ] . (10) 

The photon doistribution function of the squeezed number state |m > is descibed by the formula 


"Pn = | detCl' 1 exp [Re(d7?*C _1 d) - |d| 2 ] '■ ^ ™ m , ' • (11) 

Here m is the label of the state, whereas n is a discrete vector variable. 2Nx2N - matrix R and 
2TV- vector L are expressed now in terms of blocks of matrix ft and vector d as follows, 


R = 


( 


r 1 *? 

-r 1 


-r 1 

-v'C 1 



= R.( 



( 12 ) 


4 Even and Odd Coherent States 


The one-mode even and odd coherent states have been introduced in Ref.[7].The polymode even 
and odd coherent states have been introduced in Ref. [5] .The squeezed and correlated even and odd 
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coherent states have been introduced and studied in Ref.[4].We will discuss the photon statistics of 
the light in these states which are also called Schrodinger cat states [8].The multimode Schrodinger 
cat states are defined by the relation [5] 

| A± >= jV±(| A > ± | -A >), (13) 


where the multimode coherent state | A > is 

| A >=| ai,a 2 ,Q3, >= #(A) | 0 >, 


(14) 


and D(A) is the multimode displacement operator creating coherent state from the vacuum.The 
normalization constants are 


JV + = 


N. = 


iM 

e a 


2 yjcosh I A | 2 

|A | 3 
e 3 

2 yjsinh, | A | 2 


(15) 


where complex number A has the form 

| A | 2 =| «1 | 2 + I C* 2 | 2 + + I a n | 2= ^2 I °ri 


(16) 


m=l 


The photon distribution function has the form [5] 


P+(n) 

P-(n) 


|c*i l 2nt | &2 l 2na ... | an l 2n " 

(n x !)(n 2 !)....(n n !)cosft | A | 2 

| Ql | 2ni | l 2na - I On l 2n " 

(n l !)(n 2 !)....(n n !)sinfi | A | 2 


nx + n 2 + .... + n n = 2k, 
ni + n 2 + .... + n n = 2k + lj 


and the photon means corresponding to these distributions are 

< A + | n, | A + > = | cti | 2 tank | A | 2 , 

< A_ | ni | A_ > = | ai | 2 coth | A | 2 . 


(17) 


(18) 


The photon number dispersion matrix has the matrix elements 

a tk = I a * | 2 | ak | 2 sec ^ 2 I A | 2 + | o« | 2 tonh | A | 2 £,fc, 
a~ k = — | a, | 2 | Qfc | 2 cosech 2 | A | 2 + | a,- | 2 coth | A | 2 S ik . 

(19) 
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5 Squeezed Schrodinger Cat States 

Let us find out the photon statistics of squeezed polymode Schrodinger cat state labeled by the 
complex N -vector /?.To do that let us difine transition amplitude from the polymode squeezed 
and correlated state |/3 > to the polymode photon number state | n > 

T.W = W) JL (,-V - d)) , (20) 

where 

•W) = (det C)-» exp [|/VC ] V + ^(d* - ^rM) + - i|d| 2 - i|/?| 2 ] . (21) 

Then the photon distribution function for polymode light in the squeezed Schrodinger cat state 
(even and odd) is given by the formula [4] 

v n(0) = NiiP) [jT„(/?)| 2 + \T n (-0)\ 2 ± (T;{0)T n (-0) + T n (/3)TZ(-0))\ . (22) 

If the shift parameter d = 0 the formula is simplified 


Kifi) = 4 NlV n ((3), 


(23) 


if we have equality 

and for even states 
if one has 


N 

5Z n,- = 2 k, 

t = l 

K+M = o, 

N 

J2 n i =2k + 1, 

t = l 


(24) 

(25) 

(26) 


where < P n (/?) is given by the formula (3), For the light in the odd squeezed Schrodinger cat state 
the photon distribution is 

= WlPM, (27) 

if the indices satisfy the equality (26) and 


^2*(/?) = 0, (28) 

if the indices satisfy relation (24).Thus the squeezed Schrodinger cat states if the shift parameter 
is equal to zero have highly oscillating distribution function.The influence of shift parameter 
decreases the oscillations of the distribution function. For Hermite polynomials the following sum 
rule may be found [3] 


oo oo \ni \ns 

„7^0 n^o ni! «a! * ■ * n N \ > 

= [det (A£ X R + l2(v)] 2 exp £-z (A£ X R + I 2 a () _1 E x AzJ 
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Here z = (z t , z 2 , ...z 2J v), the 2Nx2N matrix E r is the 2iV- dimensional analog of the Pauli matrix 
cr s , and the diagonal 2Nx2N - matrix A has the matrix elements A, in j-th and (N+j)-th rows. 
Let us consider the one mode case.Then the formula for the photon distribution function in terms 
of Hermite polynomials of two variables may be expressed in terms of usual Hermite polynomials 

_ ^ (r 2 -4d)* ^ ( 4d — 1 \ k _n! 

Vn Vo (2T + 4d + l) n \>/T* -4 d) [(n - fc)!] 2 fc! 


( (T + 1)^ ~4~ [g’pp ~ <Jqq ~ 2t(7 w ] Z 


p)l' 


H k j ii + Ul + Z Z z 1 1 (30) 

\{(2T + 4d 4- 1) [o'pp — <jqq — 2t<T M ]} 

Here the parameters <rpp,<r,,,< 7 p, are matrix elements of the quadrature dispersion matrix M, d is 
determinant of this matrix and T is the trace of the matrix. The complex number z is determined 
by the relation 

z = -j=(< q> +i <p>). (31) 


Formula (30) can be used also to illustrate the generalized uncertainty relation (for the Gaussian 
states). Indeed, it is obvious that the probability to find n photons must be nonnegative. On the 
other hand, all but one terms in the right-hand side of eq. (30) are positive independently on the 
concrete values of the parameters determining the quantum state. The only exception is the term 


( 


4d — 1 
v/T 2 - Ad 



Consequently, to guarantee the positiveness of the photon distribution function for all conceivable 
combinations of the parameters one should impose the restriction d > ^ .This inequality is the 
Schrodinger uncertainty relation (see, [6], [12] ). 
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